ITEP-TH-36/99 



PHYSICS OF THE POWER CORRECTIONS IN QCD. | 



F.V. Gubarev, M.I. Polikarpov 

Institute of Theoretical and Experimental Physics, 
B. Cheremushkinskaya, 25, 117259 Moscow 
E-mail: Fedor.Gubarev@itep.ru, polykarp@vxitep.itep.ru 

0> ■ V.I. Zakharov 

On I Max-Planck Institut fiir Physik, 

Fohringer Ring 6, 80805 Miinchen, Germany 
bJO, E-mail: xxz@mppmu.mpg.de 

< 

OO ! We review the physics of the power corrections to the parton model. In the first 

part, we consider the power corrections which characterize the infrared sensitiv- 
ity of Feynman graphs when the contribution of short distances dominates. The 
^N) ' second part is devoted to the hypothetical power corrections associated with non- 

^ ' perturbative effects at small distances. 

(N ■ 

ON 

(N ■ 1 INTRODUCTION. 

00 : 

, General remarks on the power corrections. 

ON 

■ We consider QCD and processes determined by physics at short distances. Which means that there is a 

I generic large mass scale, Q ^ ^qcd where Aqcd is the position of the Landau pole in the coupling: 



For example, Q may stand for the total energy in the e+e annihilation into hadrons or the 4-momentum 
of the virtual photon in deep inelastic scattering (DIS). 
. Then one can use the perturbation theory and predictions for a physical observable O are given by a 

$—( ' perturbative series n. 

a ' / OO \ 

< O > =< O >parton model I 1 + """^(Q^)" 1 • (2) 

Now, we reserve for powers of Aqcd/Q as well: 

/ OO OO \ 

<0 > =<0 >parton model { ^ + a„as{Q^f + ^ 6„(Aqcd/Q)" (3) 



n=l n=k 



These are the power corrections. 



"Lectures given by V.I. Zakharov at the Winter School of physics of ITEP, Snegiri (Moscow Region), February 1999 
''In fact, the perturbative corrections may modify the parton model by a powers of as{Q^) as an overall factor as well. These 
are so called anomalous dimensions, best known from the example of moments from structure functions in DIS. For simplicity 
we consider the case of zero anomalous dimensions. 
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First of all, the power corrections appear to be a pure nonperturbative effect. Indeed, on one hand we 
have 



^qcd \ ^ _ 



cxpi-k/2boasiQ')). (4) 



On the other hand, the function exjp(— const/a) with a positive const is a classical example from the math 
courses of a function which has a trivial Taylor expansion at a = 0: 

exp{— const /a)\a=o = + • a + • + ... (5) 

since the function itself and all its derivatives vanish at a = 0. Thus, this function, being a non-zero, vanishes 
identically as a perturbative expansion, which is the expansion at a = 0. 
The interest in power-like corrections originates from various sources: 

(i) Nowadays, one may say that the perturbative QCD is trivially correct and, if it were possible, the 
best thing would be just to subtract it out and proceed to non-perturbative pieces. In particular, one may 
say that the physics of the confinement is encoded in the power corrections, not the perturbation theory. 

(ii) In some cases, the accuracy of theoretical fits to experimental data require for an account of the 
power corrections. For example, error bars on measured values of as are affected by the power corrections. 

{iii) More pragmatically, one could say that we are brought to consider the power corrections by the 
logic of the development in the field. 

If it is at all possible to distinguish between these motivations, we will belong rather to the first line. 
Namely, we will assume, explicitly or implicitly, that the onset of the power corrections at some so to say 
moderately large signifies new physical phenomena. At the end, we shall see whether we are in fact 
justified to assume so. 

There is another important aspect of the power corrections. While calculating the perturbative ex- 
pansions is a well defined procedure in QCD, at least as a matter of principle, the definition of the non- 
perturbative terms is close to saying that these are unknown terms, the rest of the amplitudes upon subtrac- 
tion of the perturbative contributions. In other words, working with the power corrections relies to a great 
extent on intuition and heuristic models. Nevertheless, we shall be mentioning sometimes the "standard 
picture" of the nonperturbative physics. What could this mean if we a priori know that no precise form of 
the nonperturbative fluctuations is assumed? Still, there is a content to the notion of the standard picture. 
Namely, the standard assumption is that the nonperturbative fields are soft. In other words, the typical size 
of the nonperturbative fluctuations is of order Ag^^. Later, we will challenge this picture to some extent. 

Finally, let us mention that actually working with two inflnite series as indicated in (|^) would be awfully 
difficult in practice. Thus, in reality one is always relying on a kind of a truncated series like: 

<0 > =<0 >parton model (l + ttiasiQ^) + bk{AQC D / Qf) (6) 

where k characterizes the leading power correction. The assumption behind this truncation is that the power 
corrections are somehow enhanced numerically. There is no proof of this assumption but it is an indispensable 
ingredient of any phenomenology based on the power corrections. Historically, this assumption worked very 
well in case of the so called QCD sum rules (Shifman et. al, 1979). It might well fail, however, in other cases. 
If there exist high-precision data one may try to check this assumptions varying the number of terms in the 
perturbative expansion kept explicitly and watching how the fitted value for the power corrections depends 
on this. So far, this procedure was implemented in the most careful way in (Kataev et.al, 1998) in case of 
DIS. 



Outline of the lectures. 

The power corrections are a hot subject in QCD. In particular at this School there will be another course, 
given by Yu. L. Dokshitzer also devoted to a great extent to the power corrections. To avoid overlap, this 
lectures will consist of two parts. The first one is about power corrections within the approach based on 
the operator product expansion (OPE) which goes back to QCD sum rules and is about 20 years old. The 
second part, to the contrary, describes an unconventional (and hypothetical) source of power corrections, 
that is short strings. This part is based mostly on the original work (Gubarev et. al., 1998). 
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It is worth adding that the unification of these two parts, separated by almost 20 years in terms of the 
original papers, under the same title is not artificial at all. The point is that the power corrections associated 
with the short strings hopefully solve some outstanding problems left over from the OPE-based approach. 



2 POWER CORRECTIONS AND SOFT VACUUM FIELDS. 



Correlation functions. 

It is clear that if we take the limit — > oo literally we would be left with parton model, with no corrections 
whatsoever. Moreover, the power-like corrections would die fast. Thus, our general strategy will be to start 
with large so that as(Q^) is small numerically. This is needed to have some control over theoretical 
calculation since QCD is simple only at short distances. However, then we will move down in until reach 
so to say moderate where the corrections become sizable. This is most interesting region for us. Indeed, 
we are still able to sort out various corrections since they are small compared to unity. On the other hand, 
we may hope to distinguish between various mechanisms of breaking the asymptotic freedom. In language 
of distances, wc will start with r ^ and then proceed to "moderate" r. Everywhere we understand that 
generically, r ^ '^/Q- 

For the sake of definiteness we will concentrate on correlation functions I\.j{Q^) and quark- antiquark 
potential V{r). Let us introduce these quantities in more detail. 

At first sight, it would be most natural for QCD studies to consider liadrons themselves. However, 
observing hadrons wc would not find much quarks at short distances since they are predominantly at a 
characteristic distances of order Ag^^,. And physics at such distances is governed by a large where wc 
do not have reliable theory. To ensure that quarks do not fly away one has to resort therefore to an external 
source of quarks such as electromagnetic current and consider unphysical kinematics with space-like total 



momentum of quarks q, 
for time of order 



Q ^ ^QCD- Then according to the uncertainty principle quarks can exist 
r - (7) 



which is small if is large. For consistency, after such time the quarks are to be absorbed by another 
current, see Fig. 1. 
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Figure 1: a) Correlator of currents in the parton model approximation, b) plane. 
In the field theoretical language, we are considering in fact a correlation function Iij{Q'^): 



a^e'"- (0|T{i(ar),j(0)}|0), 



(8) 



where the current j may have various quantum numbers, like spin, isospin and for simplicity we do not 
indicate these quantum numbers, i.e. suppress the Lorenz indices and so on. 

The basic theoretical ingredient is that Ii{Q^) at large Q-^ can be calculated in the parton model ap- 
proximation: 

lim IIjXQ^) = nj(Q^)parton model (9) 

On the other hand, by using dispersion relations Iij{Q'^) can be expressed in terms of the absorptive part 
which is non- vanishing only for time- like total momentum, > 0: 



(10) 



The imaginary part is directly observable, provided that the current j is a physical one. In particular, in 
case of the electromagnetic current, j — jei the imaginary part in Eq. ( p^ ) is proportional to the total cross 
section of e"'"e~-annihilation into hadrons: 

ImU, , is) = const ""^'^^'^ ^ """^^ (11) 

Upon substitution of (|Tl|), the Eq. (|lo| ) becomes a sum rule. Indeed, Iij^^{Q^) is calculable and the same as 
for free particles plus small radiative corrections, while Iml\{s) is observable. 

There is one more technical point to be mentioned. The dispersion relations (|l^) suffer in most cases 
from ultraviolet divergences which could be eliminated at a price of subtractions. But then there would 
appear an arbitrary polynomial. There is, however, another way to deal with this problem (Shifman et.al., 
1979). Namely introduce 

in the limit where both and n tend to infinity so that their ratio Af^ = /n remains finite. Then it is 
easy to see that the weight function in the dispersion relation becomes: 

while all the subtraction constants are removed by the differentiation (|l^). 
Let us give a realistic example of the sum rule in the p-meson channel: 



J Ri=i{s)exp{-s/A'P)ds w ^1 



a.(Af^) 



(14) 



where is considered to be large enough so that as{M'^)/Tr is small. Moreover i?/=i(s) is the total 
cross section of e~^e~ annihilation into hadrons with isotopic spin 1=1 in units of the standard cross section 
(T(e+e^ /i+/i"). 

It is worth noting that to apply the technique considered it suffices to ensure that the time which quarks 
exist is small indeed. Apart from imposing the condition that is large (see Eq. (|^)) there exist other 
possibilities. In particular, as far as production of heavy quarks is concerned one can consider — since 
in that case (Shifman et. al, 1976) 

r ^ l/mn ■ 

Therefore, in case of heavy quarks it is convenient to consider = and get rid of possible ultraviolet 
(UV) divergences by differentiating the dispersion relations ( [To| ) with respect to at = 0. In this way 
one comes to the sum rules (Novikov et. al., 1977): 

^#1^- + (15) 

where Rc is the contribution of the current of the charmed quarks into the ratio R{s) and the integer number 
n corresponds to the n-th derivative from ( p^ while rric is the mass of the charmed quark normalized off-mass 
shell at = 0. Moreover, An, Bn are calculable numbers: 

_ 3 2"(n + l)(n- 1)! 



47r2 (2n + 3)!! 
'(3n + 3)] TT ; 
^ + 3)] 2^4 
2 /tt 3 \ r(n + 3/2)[l-2/(3n + 6)] 4n ln2 
'V^ \2~A^) r(n + 2)[l - l/(2n + 3)] 



^ 40Fr(n + 3/2)[l-l/(3n + 3)] ^ 3 

3r(n + l)[l-l/(2n + 3)] 2 47r ^ ' 
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Note that the integral over R{s) is contributed both by resonances, Hke J/'^ and by continuum cross 
section of production of particles with open charm. The bound states are to be included since their contri- 
bution is of zero order in the small as{M'^). Indeed, the properties of the resonances are governed by as ~ 1. 
Moreover it turns out that there exist such (or n in case of heavy quarks) for which the resonances 
dominate the integral over the physical cross section while the theoretical, partonic part is still calculable 
reliably since as{M'^) is small. 

It is worth noting that direct experimental information on the correlation functions is available only in a 
limited number of cases since only very few currents, like the electromagnetic current is directly observable. 
However, the lattice simulations (Chen et. al., 1993) allow to measure ^j{x) for a much wider class of 
currents. These measurements refer to the Euclidean space directly since the lattice corresponds, of course, 
to Euclidean x. 

To summarize: there are quite a few correlation functions known in Euclidean domain either from 
combined use of dispersion relations and experimental data or from the lattice simulations. The correlation 
functions are most suitable for studies of the power corrections. 

The meaning of the OPE. Gluon Condensate. 

In this section we will describe power corrections as they arise within the operator product expansion (OPE). 
The presentation goes back to the QCD sum rules (Shifman et. al., 1979). More detailed reviews can be 
found in (Reinders et. al., 1985), (Narison, 1989). 

The central physical question which brought about the power corrections was as follows. Imagine that 
we study the sum rules ( |l^ at very large and then go down with AI^. At some low the dispersive 
integral over s would be dominated by the p-meson to such extent that it cannot be matched by the quark 
predictions. Then the question is: "Who stops the Asymptotic Freedom at some critical Af^?" 

At first sight, the answer is almost trivial: the growth of the coupling at low momenta. It would be a 
common answer. However, if one addresses the problem in concrete, there arise doubts in the validity of this 
answer. For example, it is quite a common theoretical guess that the splitting between the vector mesons 
does not depend on flavor, say: 

m{p') — m{p) « nn{%lj') — m{J /ifj), (17) 

which is true experimentally. This simple-looking observation represents, however, a serious challenge to 
the wisdom that it is the growth of the effective coupling (|l|) that stops the AF at moderate mass scales. 
Indeed, if expressed in terms of an invariant quantity, s, Eq. (|l^) implies that the J/ijj is dual to a much 
larger interval of s than the p because the c-quark is heavy. In other words, the AF is violated in the p 
channel much later than in the charmonium channel, if we start from very large downwards. A direct 
numerical analysis of the sum rules (|4|), (|5|) confirms this expectation. However, the coupling should run 
as a function of an invariant quantity and is flavor blind. 

More generally, one feels that knowing the perturbation theory alone is not enough to plunge into 
the resonance physics. Therefore, let us assume for the moment that it is soft nonperturbative fields in 
the vacuum which are eventually responsible for the confinement. At first sight, having this idea does 
not help much since very little is known on the precise nature of these nonperturbative fields. The way 
out of this difficulty will be not to try to calculate the nonperturbative fields but describe them instead 
phenomenologically in terms of a few parameters. 

One can understand the trick by considering the same Feynman graphs for the Yij{Q^). Let us begin 
with the simplest graph of Fig. 2. As was explained above (see Eq. (Q)), by taking large we ensure 
that the correlation function is determined by short distances. Let us now estimate to which accuracy this 
assertion is correct. Consider to this end in more detail the route of the large Q, which is brought in by the 
current. The typical case is when the both quark lines in Fig. 1 carry momentum of order Q. "Typical" 
means favored by the phase space. But it is not forbidden that almost the whole momentum is carried by one 
line while the other quark line is soft. Then we get a problem since the soft line cannot be reliably described 
by the perturbative propagator because at large distances the perturbation theory is modified strongly. 

Instead of making guesses, how the quark propagator looks like in the infrared (IR) we reserve for an 
unknown number for this propagator. The central point is of course that it is not an unknown function but 
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rather a number. To see that it is indeed true, prepare the graph as is shown in Fig. 2. Here the black 
k 




Q 



Figure 2: If one of the lines is soft, one substitutes instead the perturbative propagator the corresponding vacuum matrix 
elements of the fcrmionic operator. 



blob in the right-hand side denotes the hard quark line and currents. Along this line, since it is hard, all 
the distances are small, x ^ 1/Q and everything is calculable. On the other hand, the soft line can travel 
distances of order 1/k ^ ^qcd which are large on the scale of 1/Q- However, it is constrained to begin and 
end up at the same point. We can substitute therefore the soft line by an unknown matrix element. This 
matrix element is proportional to the so called quark condensate, < OjgqlO > 7^ 0, which is famous for the 
fact that its non- vanishing value signals the spontaneous breaking of the chiral symmetry. If the current j is 
constructed on the light quarks, then it is easy to see that in fact the quark condensate is multiplied by the 
quark mass and the contribution of the IR region is additionally suppressed. Namely it is of the order 

STljjQ') m, < 0\qq\0 > 

n,(g2) ^ • 

Notice the factor which we get on pure dimensional grounds. This suppression is the price we pay for 
enforcing one of the lines to be soft. 

The graph in Fig. 2 is very important, however, if one of the quarks is heavy. Then we replace rrig in 
Eq. (|l8|) by mn while q still stands for a light quark. This is our first example of a nontrivial interplay of 
perturbative and nonperturbative calculations. Indeed, on one hand, we consider perturbative graph. The 
perturbation theory by itself respects the symmetries of the Lagrangian and < 0|qg|0 >pert^ and is 
negligible. However, in the general treatment of the infrared sensitive contributions which we pursue now 
the quark condensate < > is just a phenomenological number which is in fact does not go to zero for 
ruq 0. 

Coming back to consider currents constructed on the light quarks, the contribution ( p^ is small because 
mq is small and we have to go to the next order in perturbation theory to uncover IR sensitivity of the 
Feynman graphs at large Q. Thus, turn to the graph with a gluon exchange. Assuming that the momentum 
flowing through the gluon line is in fact of order Aqcd we prepare the graph as indicated in Fig. 3. Again, 




Figure 3: Treatment of one-gluon exchange graphs according to the OPE in case when the gluon line is soft. The black blob 
corresponds to the hard quark lines and the currents. All the points here are close to each other. Instead of the soft gluon line, 
denoted by a dashed line, one substitutes a phenomenological matrix element. 
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there is no reason any longer to use the perturbative expression for the giuon hne since it is soft and modified 
strongly by the confinement. The quark lines propagating short distances become a receiver of long wave 
gluon fields in the vacuum. The receiver is well understood because of the asymptotic freedom. The intensity 
of the gluon fields is measured. It is characterized by the so called gluon condensate: 

(0K(G;J2|0) ^ (OK ((i/'^f - |0> ^ (19) 

where i?", E'^ are color magnetic and electric fields and G^^, is the gluonic field strength tensor. 

The question may arise, how do we know that it is the matrix element that arises in a straightforward 
evaluation of the IR sensitive part of the one-gluon exchange graph. The answer to this question is not 
difficult: this is the simplest expression compatible with the Lorentz and gauge invariance. In other words, 
the (G^y)^ is the simplest operator which is Lorentz and color singlet. It has dimension d=4 implying that 
the IR sensitive contributions are suppressed at large as Q"**. 

What is described here in words, in its generality has an adequate formulation in terms of the Wilson 
operator expansion (Wilson, 1969). A systematic treatment of the power corrections within the Wilson 
OPE was given in (Shifman et. al., 1979). In particular, it allows to evaluate the contribution of the gluon 
condensate (|l9|) to correlation functions corresponding to various currents j. 

One may still wonder, how it could be possible to go beyond the perturbation theory through a simple 
preparation of perturbative graphs. To answer this question, let us try to evaluate the gluon condensate in 
perturbation theory. To lowest order it is given by a simple one-loop graph: 

{0\asiG;,r\0) « / d'k^, (20) 

where we did not cancel immediately in the denominator and numerator of the integrand to emphasize 
that the factor k~'^ is the gluon propagator while the factor k'^ upstairs is associated with the derivatives in 
the vertex. Also, we assumed that there are eight color gluons, as in the realistic case. From the first glance 
at the expression (20) it is clear that it diverges in the ultraviolet as A^y. However, in all the applications 



we assume instead that 

where the constant cg is to be found from the fitting procedure. And there is a rational behind this 



\asiGlS-\0) ^ cgA%cd (21) 



madness. The point is that the UV divergent part of the integral (21) corresponds the ordinary perturbative 
contribution. (Of course, there is no real UV divergence in the original Feynman graph (see Fig. 3) since 
at A; ~ Q the account for other lines in the graph provides with a kind of form factor so that the final 
result is finite.) Thus, we include only the IR contribution into the gluon condensate. Generally speaking, 
we should have subtracted then the perturbative IR contribution, included now into the matrix element of 
(GJ^j,)^ from the ordinary perturbation theory. This is to avoid a double counting, as was emphasized by 
David (1984) and Mueller (1985). There are no difficulties of principle to deal with the problem (see, e.g., 
Shifman et al (1979), Novikov et.al. (1985)). However the whole approach is becoming most interesting from 
the phenomenological point of view if the nonperturbative contribution is enhanced numerically. And this is 
the assumption which is commonly made at this point. Thus, the OPE allows to introduce in a consistent 
fashion the notion of the enhancement of the IR sensitive contributions to Feynman graphs. 

Coming back to the problem of isolating soft lines in perturbative graphs, the same one-gluon exchange 
graph can be characterized by another route of the large Q flow, see Fig. 4. Here we have two soft quark 
lines which means in turn that the corresponding correction to the right-hand side of the sum rules ( p^ is 
of the type: 

< 0\qOqqOq\0 > 

where O is constructed on the Dirac matrices 7^ in the spinor space and on the Gell-Mann matrices A" 
in the color space. Usually one assumes factorization of the 4-quark matrix element reducing it to the 
quark condensate mentioned above. This time, however, the condensate is not suppressed by power of the 
(practically vanishing) light quark masses. This is another example how the nonperturbative phenomenon 
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Figure 4: Space-time picture corresponding to the Wilson operator product expansion in case that there are two soft quark 
lines. 



of the spontaneous breaking of the chiral symmetry finds its way into the OPE based on analysis of the 
perturbative graphs. 

With account of the quark and gluon condensates the sum rules in the p-channel become: 

Rpis) exp{^s/M^)ds « x (1 + + (23) 



2 ' TT 

TT (0|a.(G°J^|0) 32 ■ 14^3 \{o\al/^qq\o)\2 
3 M4 81 A/6 



■) 



where the ellipses stand for perturbative and power corrections of higher order. In case of the charmonium 
sum rules (see (p7[)) there also appears a term proportional to the gluon condensate. 



Eqs. (23) are the QCD sum rules (Shifman et. al., 1979). The first check was to see whether the 
gluon condensate explains the difference in the duality intervals in the p- and J/ip- channels (see above). It 
does explain this difference naturally. Also in many other cases, the QCD sum rules turned to be a very 
straightforward and successful tool for orientation in the hadronic world. For our review it is most relevant 
that the sum rules provide with ample confirmation of the idea that it is the soft vacuum fields that are 
responsible for the breaking of the asymptotic freedom at moderate and signal in this way the formation 
of the resonances. One can characterize the power corrections by M^^^^ which is defined as the value of 
A/^ where the power corrections become, say, 10% of the partonic contribution. Then the typical M^^^^ 
associated with the quark and gluon condensates is 

Ml,, « 0.6 GeV^ (24) 



and its dependence on the channel for currents constructed on light quarks is not strong. The meaning of 
the A/^^jf is that the power corrections blow up at smaller A/^. Clearly, the value ( ^ ) resembles the p- meson 
mass squared. 

To summarize: the use of OPE allows to parameterize infrared sensitive contributions to the Feynman 
graphs in a systematic way in terms of vacuum expectation values. While the coefficients in front of the 
matrix elements are calculable as an expansion in as(Q^), the matrix elements encode the information on 
nonperturbative phenomena. 



The heavy quark potential. 

Another useful quantity to consider in connection with the power corrections is the heavy quark potential 
at short distances, for a review and further references see (Akhoury and Zakharov, 1998). 

Let us first notice that although the static potential between heavy quarks is obviously a fundamental 
quantity, its definition in non-Abelian theories is not straightforward at all. The point is that we are interested 
in case when the quarks are in a color-singlet state. On the other hand, the gauge invariance allows to rotate 
quarks in the color space locally and, at first sight, it is just impossible to single out a color singlet state. 
The way out is to consider the so called Wilson loop, that is the average of an operator constructed on the 
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gauge field potential A^j^ along a closed path C. Moreover, it is convenient to choose the contour C as a 
stretched rectangle C ^ r x T with large T. Then the potential V{r) is: 



V{r) = - Jim -^\ii{Tr P exp ( ig j> dTf^A^T" )). (25) 



T-*oo iT 



This definition is gauge invariant and refers to the Euclidean space. The main source of knowledge of V{r) 
are the lattice simulations (see Bali et al, (1995); Bali, (1999)). 

The most famous property of the potential V{r) is that it grows linearly at large distances r: 

lim Vir) - floor. (26) 

r — *oo 

Note that here we imply that the dynamical light quarks are not taken into account and consider the 
interaction of infinitely heavy external quarks embedded into the vacuum of pure gluodynamics. 

We are mainly interested in the potential V{r) at short distances. Let us begin our analysis with a trivial 
Yukawa attractive potential and expand it in powers of (Ar) at small r: 

VYu{r) = --e ^' = -- + CA 7^ + ^^ + - > 27 

r r 2 6 

where C is a positive constant. This simple equation carries an important message, namely, that the physics 
behind the odd and even powers of the mass A in the expansion (^^ is different. Indeed, naively we should 
have had only powers of A^ since A^ is the only mass parameter entering the Lagrangian of a massive boson. 
Thus, our idea could be that we start with a Coulomb like potential at small r and develop a perturbative 
expansion in X^r^. Eq. ( p7| ) implies that such an approach would fail because of infrared divergences. This 
can be readily checked explicitly of course. But it suffices to notice that the odd powers in ( |27| ) could arise 
only from infrared divergences which are cut off at distances of order 1/A. 

Appearance of terms non-analytical in A^ in ( [2^ ) allows to make important conclusions about the heavy 
quark potential in QCD. Generally, if a physical observable is sensitive to large distances at level of a certain 
power correction then it cannot be evaluated to such accuracy because the effective coupling is large. Rather, 
one should reserve for the strength of the correction as a phenomenological parameter. In our case, we can 
conclude from Eq. ( p7[ ) that the quark-anti-quark potential in QCD can be parameterized as: 

V{r) = ^^ + CoAQCD + C2r^AlcD + - (28) 

where C_i is calculable perturbatively as an expansion in as{r) while Co,2 account phenomenologically for 
the contribution of large distances. Moreover, Cq can be included into definition of the (heavy) quark masses 
and the message is that the power corrections to the potential at short distances start with terms (Balitsky 
(1986). Dosch and Simonov(1988)). A salient feature of Eq. (ESh is absence of a linear in r term at short 



distances. The proof was in two steps. First, we learn from (27) that the linear correction to the potential is 
not contaminated by large-distance effects and, then, conclude that these corrections vanish since the only 
parameter of dimension d=2 is the gluon mass squared and it vanishes in QCD. Thus, presence of a linear 
term at short distances in the quark potential would signal new non-perturbative physics at short distances. 



The crucial assumption behind Eq. (28) is that typical size of nonperturbative fluctuations is of order 
^QCD- This connection can readily be visualized in the abelian case. Namely, begin with an identity for the 
potential: 

V{r) = j- y'dV'Ei(r')-E2(r-fr'), (29) 

where Ei, E2 are electric fields associated with each charge. Eq. (|9|) is convenient to account for the effect 
of the change in the electric fields at large distances. Consider an example of two charges of opposite signs 
in a cavity of size R. The electric field of the charges in empty space is that of a dipole at large distances, 
E'^ ~ ar^/r ^. Now, because of the cavity the field is changed at r ~ i? and the corresponding feedback to 
the potential V{r) is of order: 

2 f°° <Pr' ar'^ 



'5^(0 - / (^7)6 - 1^ (30) 
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which is in agreement with the correction to the static potential discussed above. 

After, hopefully, explaining thoroughly that the correction to the potential V{r) at short distances is of 
order in the standard picture we are going to introduce a new notion, that is retardation effects (Casimir 
and Polder, 1948). Discussion of this point is a kind of deviation from the course to our immediate goals. 
We are in fact interested in the static potential at short distances. And then there are no retardation effects, 
of course. However, historically atom-like systems of heavy quarks were discussed at most and in this case 
the retardation effects may be important. This is the only reason to include discussion of the retardation 
effects into this review. 

Thus, Casimir and Polder did discuss the interaction of atoms with nonperturbative long-wave fields 
in QED about 25 years before the advent of QCD. What kind of nonperturbative fields did they consider? 
At first sight at least, there no nonperturbative fields in QED. They considered atom in a cavity. Then 
zero-point fluctuations are modified by the presence of the cavity at frequencies of order uj ^ R^^ where R 
is the cavity size. If in the QCD case, we substitute the cavity of size i? by a nonperturbative field of the 
size ^Qcu, the problems turn similar! 

Of course, here we present only a very sketchy view of the paper by Casimir and Polder. One starts 
with the dipole interaction (see also above): 

H,nt = -ed-E. (31) 
The shift in energy levels of atoms arises then in second order in this interaction: 

^ ''^VnkiEn — Ek + LUchar)~^Vkn, (32) 
k 



where the characteristic frequency of the nonperturbative fields is of order ujchar ^ 1/R- The result for the 
energy shift SEn depends crucially on the relative magnitude of ujchar and £"„ ~ ma^ ^ a/as where as is 
the Bohr radius. Namely, if uJchar 3> En, then 



a2 

5Eo - ad^E^i? - a-S-, R<s:aB/a (33) 



since d^ ~ a^, and E,^ ^ R 



The shift (33) corresponds to our estimate given above in the language of the classical electrodynamics. 
In this case, the potential picture applies to the evaluation of the energy shifts. 

On the other hand, if i? ^ as /a then cjchar in the energy denominator in ( |3^ ) can be neglected and 

3 

6Eo ^ ad^E^^ ~ ^, i?»aB/a, (34) 

where we used the estimates of E^,d^ given above and substituted En ~ a/as- 

The estimate ( |3^ ) is in clear violation of the potential picture. Moreover, Eq. (^ ) could be interpreted 
by saying that if the distance between the particles is much smaller than aR then the potential is 
replaced by a potential ^ Note also that the emergence of the scale R ^ UB/a in the equations above 
can be understood as an effect of retardation. Indeed, the time needed to communicate with the distances 
of order of the cavity size R can be called the retardation time. Tret ^ R- For the potential picture to be 
valid. Tret is to be smaller than the revolution time which is of order, Trev ^ clb/v ~ aB/ot. The potential 
picture becomes distorted once Tret ~ Trev 

In QCD case one considers (see Voloshin (1979), Leutwyler (1981)) atom-like systems of heavy quarks 
Q such that 

as{MH)-MH > ^QCD (35) 



"^But this is true only as far as rough estimates are concerned. Rigorously speaking, there is no potential whatsoever corre- 
sponding to the shifts obtained in this way. This was emphasized much later, in connection with QCD, by Voloshin (1979) and 
Leutwyler (1981). 
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where Mh is the mass of the heavy quark and the condition ( ^5| ) is that the Bohr radius is much smaUer 
than ^QQjy- Again, the dipole interaction is relevant to describe interaction with soft gluon fields: 

H^nt - - VSr(t? - t«)d • E^ (36) 

where tf {i — 1, 2) are generators of the color group and refer to the quark Q and anti-quark Q in the 
quarkonium while denotes the soft gluonic field in the vacuum. 

Moreover, one assumes that uJchar ~ ^qcd- intensity of the gluonic fields, it is characterized 

again by the gluon condensate < as(G°y)^ >. The resulting shift of the energy levels of the S'-states depends 
strongly on the principal quantum number n: 

Mh [MHCFasr 

where the numbers e„i depend on the quantum numbers of the level, n, Z, and are of order unity, e.g., 
£10 ~ 1.5. The growth with n is due to the growth of the size of the atom-like state. However, to apply ( |37| ) 
one should assume that the retardation effects are already very important. For a concrete value of the heavy 
mass Mh the balance between these requirements can be quite delicate and we refer the reader to a review 
by Yndurain (1998) for details and references. 

To summarize: according to the standard QCD there is no linear correction to the Coulomb potential 
at short distances. This is in variance with naive unification of the linear potential (dominating at large 
distances) and Coulomb-like potential (dominating at small distances). 

Moreover, because of the retardation effects the effect of the leading correction in V{r) would be 
washed out in heavy quarkonium and the nonperturbative corrections to the Coulombic potential at short 
distances would start with even smaller terms of order . 

Other techniques: infrared renormalons, infinitesimal gluon mass. 

Notice that while discussing V{r) we did not resort to the OPE. It is not by chance of course. The point 
is that no matter how simple the notion of the potential V{r) could appear, the OPE cannot be applied 
directly to evaluate V{r). The reason is that we consider the static potential which means, potential energy 
averaged over large period of time. Thus, although we consider small spatial distances r, the potential V{r) 
is sensitive to large intervals in the time directions. 

It is not the only case when the OPE does not apply of course. Moreover, as was mentioned above, 
the OPE applies actually in the Euclidean space and only in some cases, when the analytical properties are 
simple, the results can be continued to the Minkowski space. 

Thus, it is desirable to develop techniques to evaluate the power corrections directly in the Minkowski 
space. And there exist such techniques. In particular, we will discuss here infrared renormalons and in- 
finitesimal gluon mass. We do not have time to cover the topics in detail. Instead we will illustrate the 
technique by two examples, rederiving the results which we already know. A more detailed exposition and 
further references can be found in reviews, see, e.g., Dokshitzer et. al. (1996), Akhoury and Zakharov (1997), 
Zakharov (1998), Beneke (1998). 

First, we will consider the power corrections to the potential V{r) by using the IR renormalon technique 
(Aglietti and Ligeti, (1995); Akhoury and Zakharov, (1998 )). Consider to this end the one gluon exchange 
potential 

N ^ f , exp(ik • r) 

= (^47ra.(k^) \ (38) 

where Cp is related to the color indices and Cp — 4/3 in the realistic case of the S'[/(3) color group. Note 
also that running of the as(k^) incorporates the leading logarithmic corrections. To find the renormalon 
contribution we rewrite (|l|) identically as 

^ \ — 2crfco 
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Next, substitute (^9|) into ( Pq ) and perform integration over directions of k to get 

2 p sin(fcr) ^ ^ 



Moreover, since we are interested in V{r) at small r we expand in kr: 

V{r) = daA^^^}j J dk{k~^"''° + fc2-2-fco^2 ^ ^ ^4^^ 

Now, integrating over k introduces what is called renormalon poles at 

1 3 

^poie = 2^0, ■■■ (42) 

The next integration, that is over cr, becomes undefined as a result of these poles. One can resort, say, 
to defining the integrals as their principal values. But this prescription is openly arbitrary and one is to 
reserve for an overall rescaling of the pole contributions. As a result, we get as corrections to the Coulombic 
potential at short distances: 

5V{r) = coKQCD^C2Kl^cDr^ + - (43) 
reproducing the result speculated upon above. Since we expand in kr and k ^ ^qcd^ ^'i- (El) valid at 

Since we were able to determine the leading powers for the power corrections by means of the renormalons, 
one could think that the OPE is no better than the renormalons. It is, unfortunately, not true. Namely, 
renormalons do not allow, generally speaking, for a systematic expansion of the coefficients in front of the 
power corrections in small coupling, like as(k^). Generally speaking, all orders of perturhative expansion 
collapse to the same order contribution if projected onto the power corrections. 

Let us illustrate the point by the same example of the quark potential. In fact three first orders of the 
perturbative expansion in front of 1/r were explicitly found, see Peter (1997): 

V{k^) ^ - "^^P^^MS , ^1^. ^^(2.583 - 0.278n/) + ^^(39.650 -4.147n/ + 0.077n^)j (44) 
Moreover, in the approximation of the one-loop /3-function, 

a^ik') = -^AgcD-r^a.(fc^). (45) 

ZOo dAqcD 

By differentiating n times as with respect to In Aqcd we immediately find contribution of (n + 1) renormalon 
chains, associated with the a"^^(fc2) in the perturbative expansion. 

Thus, we have a unique possibility to compare the contributions of one-, two- and three renormalons: 

5V{r) « (co?isi)A|c7^r2(l + 1.1 + 6.0 + ..). (46) 

We see that there is no convergence in sight, even numerically ^ 
The general rule is that if perturbatively 

Vik") = {A'Kaje)f{as) 

then 

5V{r) = vokQCDf^'ias - l/2&o) + 1^2 A^^^r^ /""(«. = 3/2&o) (47) 



''strictly speaking, we should have kept further terms in the /3-function since we consider a multi-loop effect. Thus, the example 
is still rather illustrative than rigorous. 

12 



where wo,2 are constants and / is the Borel-miproved expansion, i.e. the series with the expansion coefficients 
= a„/n!. 

Another technique which can be used directly in the Minkowski space is the introduction of a (fictitious) 
gluon mass A, A ^ and tracing terms non-analytical in A^. In fact, we used this technique when analyzed 
the Yukawa potential. Now we would like to emphasize its generality. Moreover, the idea is similar in fact 
to that underlying introduction of the gluon condensate. 

Indeed, let us consider again the gluon condensate perturbatively (see Eq.(^0|)) but this time in case of 
a finite gluon mass (see Fig. 5): 

It diverges wildly in the ultraviolet, as discussed. Let us define the gluon condensate, however, as the 
non-analytical in A^ part of the perturbative answer (Chetyrkin and Spiridonov, 1988): 

(0|«.(G;5J^|0) =^/°°(^^f^--^A^lnA^. (49) 

Note that there is no difficulty in practice to pick up the non-analytical in A^ terms. To this end, it is 
sufficient to differentiate (^) twice with respect to A^. 

Moreover, to finally get rid of the gluon mass (which is not a pleasant sight for a theorist's eye) we make 
a replacement: 

a^A^'lnA^ ^ CiAgcD (50) 

where C4 is an unknown coefficient ^ The central point is that if we evaluate various correlation functions 
(^ and isolate the terms A'^lnA^ we would reproduce the sum rules (23). Indeed, as we explained above 
the gluon condensate in the sum rules (23) parameterizes the infrared sensitive part of the Feynman graphs 
associated with a soft gluon line. Picking up terms non-analytical in A^ is the same good for this purpose 
since the non-analyticity in A^ can obviously arise only from soft gluons, fc ~ A. 

So far we have not got any new result, though. The real advantage of introducing A 7^ is that 
calculations can be performed now in Minkowski space as well and apply for this reason to a much wider 
class of observables than the OPE underlying the QCD sum r ules (H) 0. The link to QCD is again through 
the bald replacement of the non- analytical in A^ terms by corresponding powers of Aqcd- 

asVX^ -> ciAqcd 
asX^lnX^ ^ C2AlcD ■•■ (51) 

where Ci_2 are some coefficients treated as phenomenological parameters. The gluon condensate appears now 
merely as one of the terms in this sequence. 

The phenomenology based on such rules turns successful (see, e.g., Webber (1999) and references therein). 
The problem is not so much a lack of success but rather too much of overlap (Akhoury and Zakharov, 1996) 
with old-fashioned hadronization models, like the tube model. 

To summarize: there exist simple techniques which allow to parameterize the infrared sensitive parts 
of the Feynman graphs directly for observables in the Minkowski space. The success of this, most naive 
approach reveals that at least in the cases when the technique applies the nonperturbative effects reduce to 
a simple amplification of infrared sensitive contributions to the Feynman graphs. 

" The technique with A 7^ is in fact close to the technique based on the infrared renormalons, mentioned above.. In the 
context of the gluon condensate, the infrared renormalons were considered first in (David, 1984) and (Mueller, 1985). 
■'^There is a price to pay, however. Infrared renormalons if applied directly in the Minkowski space do not allow for model 
independent relations between power corrections to different observables. The reason is that the coupling as refers now to 
infrared region and is, therefore, of order unity. As a result, all orders in as are equally important, see (Zakharov, 1998). This 
is the same "collapse" of the perturbative expansion on the level of the power corrections illustrated above on the example of 
the potential V{r). 
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Direct instantons 



This simple picture is not universally true, however. First examples when it fails were found as early as in 
1981 (Novikov et.al.) Namely, there exist channels where the infrared sensitive corrections described above 
cannot be the whole story. To give counter-examples, one can either rely on the analysis of experimental data 
or on theoretical methods which would allow to evaluate the power corrections without using the OPE. Both 
ways were exploited to demonstrate that there exist corrections which can be very important numerically 
and which go beyond the picture described above. 

In particular, in the 0^-gluonium channel one can establish a low energy theorem relating the value of 
the correlation function at = to the same gluon condensate < q;s(G°j^) >. Moreover, one can utilize 
this knowledge to subtract the dispersion relations and convert the 11(0) into a power correction at large 
momenta. As a result the following sum rules arise: 

J ImG{s) cxp {-s/M^)ds/s « 

(0|«,(G°J^|0) / 27T^ 16^^ W 

where 

G(g2) ^ijd^x (0| T{{Gl^{x))\ {GlMf} |0). (53) 

To reiterate, the power correction proportional to the gluon condensate originates here from two sources. 
First, there is a standard OPE correction (see (|23|)) and, second, the one evaluated via a low energy theorem 
specific for this particular channel. The correction which is not caught by the standard OPE is about (20-30) 
times larger! 

Thus, if we characterize the scale where the asymptotic freedom gets violated by the power correction 
by the value M^^^j as is mentioned above, then M'^^^^ differs drastically in some channels: 

(Af2^,,)p_,„ « 0.6 GeV^ (M2„Jo+ gi^omu^ « 15 GeV^ (54) 

We see that the proof of the low-energy theorem brought a proof of existence of qualitatively different scales 
in the hadron physics (Novikov et.al., 1981). 

One may call the 0+-gluonium channel exceptional because of the appearance of this new correction. 
In case of the 7r-meson channel, it was also possible to show that the standard OPE corrections are not big 
enough to match the contribution of the pion into the dispersive part of the sum rules. Moreover, there 
arises a kind of hierarchy of exceptional channels, or of M^^^^. 

This hierarchy could be explained qualitatively in terms of transitions of the corresponding currents 
directly to instantons, see Fig. 5 (Novikov et. al., (1981), Geshkenbein and loffe (1980)). At large Q the 



G{M 



2\ 

J parton model 




Figure 5: Direct instantons. One substitutes instanton fields, both bosonic and fcrmionic into the currents and integrates over 
the instanton sizes. 



direct instanton contribution dies off fast, like (Q^) '^■^ . However, at intermediate these corrections may 
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become important first because of a big overall coefficient. Unfortunately, one cannot go far with clarification 
of the situation by using only analytical means. Nowadays one relies on the model of instanton liquid which 
allows for a much more quantitative treatment of the instanton effects (for a review and further references 
see (Shuryak and Schafer, 1997). The model appears to be successful phenomenologically. 



3 BEYOND THE OPE. 
Elusive effects of conflnement 

Looking backward, it still remains a mystery whether any specific confinement effects are revealed through 
the power corrections discussed so far. Indeed, consider the vacuum of pure gluodynamics. It is known 
from lattice measurements that external heavy quarks are confined by this medium. On the other hand, the 
effects included into the sum rules so far do not seem to encode the confinement. Indeed, the perturbative 
QCD resembles ordinary bremsstrahlung in QED. The gluon condensate, as well as other newly found power 
corrections can be detected by introducing a fictitious gluon mass which is not sensitive to the non Abelian 
nature of gluons at all. Finally, instantons are known not to ensure the confinement either (Chen et.al., 
1998). 

In an attempts to find power corrections related more directly to the physics of confinement one can 
turn to the Abelian Higgs model (AHM) which underlies the dual superconductor model of the confinement 
(Mandelstam (1974), Nambu, (1974); Polyakov (1975), 't Hooft (1976)). The basic idea behind the model is 
that the properties of the QCD vacuum are similar to the properties of an ordinary superconductor. Indeed, 
if a pair of magnetic charges is introduced into a superconductor, the potential energy of the pair would 
grow linearly with the distance r at large distances: 

lim V{r) « CToo • r (55) 

r — *oc 

where (Too is the tension of the Abrikosov-Nielsen-Olesen string. In QCD, a similar phenomenon was postu- 
lated to happen, with a change of magnetic charges to (color) electric, or dual charges. 

If one introduces a pair of external magnetic charges into the vacuum of the AHM model in the Higgs 



phase then the potential grows at large distances, see Eq. (55). The scale of distances is set up by the inverse 
masses of the vector and scalar fields, rriy^g. This growth of the potential is due to the Abrikosov-Nielsen- 
Olesen strings. 

Consider now short distances, r <C Wys- Then the Coulomb like interaction dominates. However, there 
is a stringy correction to the potential at any small distances (Gubarev et. al., 1998): 

lim V(r) = - ^ + aor (56) 

where Qm is the magnetic charge. The ANO string is a bulky object on this scale and is not responsible for 
the linear correction. Instead, the stringy potential at short distances is due to an infinitely thin topological 
string which connects the magnetic charges and which is defined through vanishing of the scalar field along the 
string. Thus, it turns out that at least in this model the confined charges learn about confinement already 
at small distances because of the short strings which are in fact seeds for future confining ANO strings. 
Amusingly enough, it demonstrates that at short distances a dimension two quantity is not necessarily the 



gluon mass squared but could be a string tension as well. Eq. (56) is the main result which we are going to 
substantiate in this part of the lectures. However, we need to make a few steps before we can explain (|5^). 
Also, we shall spend some time later to discuss the validity of the OPE. 



Field theoretical solenoid. 

A key element in the conjectured mechanism of the quark confinement is the Abrikosov-Nielsen-Olesen string 
(Abrikosov. 1957), (Nielsen and Olesen 1973). Although eventually it will not be this string that captures 
our attention, it is useful to begin with a brief review of the ANO string. 
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The strings are classical solutions to the Abelian Higgs model (AHM). The model describes a gauge field 
interacting with a charged scalar field $ as well as self-interactions of the scalar field. The corresponding 
action is: 

S = J d^x^^^F^^, + ^\{^-^A)<^>\^ + \xm'~v'r] (57) 

where e is the electric charge, A, j] are constants and F^i, is the electromagnetic field-strength tensor, _F)j^ = 
d^Ay — d^Afj,. The scalar field condenses in the vacuum, ($) = rj, and the physical vector and scalar particles 
are massive, my = e^rf , mj^ = 2Xri^. In the perturbative regime, interaction of the massive scalar and vector 
particles can be readily calculated order by order. 

There exists, however, a topologically non-trivial stringy solution to the classical equations of motion 
which possesses a cylindrical symmetry and is characterized by a finite energy per unit of its length. In 
a way, it realizes an analog of a solenoid in field theory. The key element is to look for a solution with a 
non- vanishing electric current 

^ e(<i>*(a^<i>)-(a^$*)<i>). (58) 
Moreover, if one chooses the scalar field of the form, 

$(r) = e'-^fip) (59) 



then the electric current is circular and reminds, therefore, the solenoid current. The value of n in Eq. ( ^91) is 
integer for the field $ to be a unique function of the coordinates. Since the current has only a nonvanishing 
j0 component, the functional form of the matching vector potential is 

A(r) = e^^- (60) 

The central question is whether it is possible to ensure finiteness of the energy (per unit length) with the 
ansatz (|59|), (|6C|). The behavior of the fields at p ^ and p ^ cx) is crucial at this point and it is easy to 
check that the conditions 

/(O) - A{0) (61) 

/(oo) = V 

A(oo) = - 
e 

suffice to make the energy finite. 

From these conditions alone, one can derive the magnetic flux carried by the ANO string. Indeed, 

HdS = d) A^dx^, = — n, (62) 

where to evaluate the latter integral we used the asymptotic value of A^, see Eq. (^l]). Thus, the flux is 
an integer of a minimal magnetic monopole charge. It is obvious then that the strings can end up with 
monopoles. 

Explicit form of the functions /(p), A{{p) introduced above can be found by solving the classical equations 
of motion: 

A,, = ^ ($*9^$ - $9^$*) - e'A^ (63) 

V^* + 2ieAf,^^'^ - e^Af.A'' = Xif^ - A$|$p 

In particular, one concludes from these equations that the function f{p) approaches its asymptotic value at 
p ^ CO exponentially, as exp(— m^rp). The magnetic field falls off exponentially, as exp( — myp). Notice, 
however, that the vector potential A^ falls off only as an inverse power of p. There is no field strength tensor 
associated with the asymptotic value of A^. As for the full functions f{p),A{p) they can be obtained by 
numerical methods. 



16 



Details of the solution depend on the ratio of the masses, my and mn- Two limiting cases are of special 
interest. In the so called Bogomolny limit, 



mv = ruH = rn, (64) 

the equations simplify greatly and can be solved in fact analytically, as a series in mp. Another useful case 
is the so called London limit, 

m'jj > my (65) 

In this limit, the scalar field is frozen as = 77 and one can neglect in most cases the fluctuations of the field 
$ around its vacuum value. However, the string tension, or energy per unit length becomes logarithmically 
divergent in this limit: 

' (66) 



niy 

Indeed, without the Higgs field one cannot construct a solution with a finite string tension. 

From the physical point of view, the most important manifestation of the ANO string is the linearly 
growing potential energy for an external monopole-antimonopole pair (|5^), when the distance r is much 
larger than the inverse masses, m^"^^. Thus, the magnetic charges are confined. This property of the 
Abelian Higgs model is the basis for the dual superconductor mechanism of the quark confinement. Namely, 
one assumes that the QCD vacuum is similar to the vacuum of the AHM in the sense that a scalar field with 
non-vanishing magnetic charge condenses in the QCD vacuum. Then the (color) electric charge is confined. 

Note that once we allowed for point-like monopoles, we should have considered the energy of the Dirac 
string as well. Let us discuss this issue on another example, that is compact U{1) gauge theory. 

Compact Photodynamics. 

In this section we will outline the paper (Polyakov, 1975) which, from the point of view relevant to the 
present lectures, demonstrates that there can exist a very different source of nonperturbative effects than 
those we discussed so far. 

The action of the theory we are going to consider is very simple: 



^2 I d'^n^ (67) 



where is the abelian field strength tensor, i^^^ = d^^Ai, — d^A^. The action (|67| ) is that of free photons 
and at first sight nothing interesting can come out from this theory. In particular, if we introduce external 
static electric charges as a probe, their potential energy would be given by one-photon exchange without any 
corrections. 

However, we shall see in a moment that, in a particular formulation, the theory admits also magnetic 
monopoles. Hence, a few preliminary words on the monopoles. Monopoles have magnetic field similar to the 
electric field of a charge: 

H = (68) 
Then the flux of the magnetic field through a surface surrounding the monopole is: 

$ = <ZM (69) 

On the other hand, because of the equation div H = the magnetic fiux is conserved. Thus, the magnetic 
monopole cannot exist by itself and one assumes that there is a string which is connected to the monopole 
and which brings in the flux. Moreover, to make the string invisible one assumes that the string is infinitely 
thin. Finally, to avoid the Bohm-Aharonov effect one imposes the Dirac quantization condition, 

e (j) Afj,dxi^ — e J H. ■ ds ~ 2tt ■ n , or qm — — n (70) 
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Also, we ask for the energy (or action) associated with the Dirac string to vanish. Only then energies of the 
electric and magnetic charges are similar. We shall return to discuss the issue of the energy of the Dirac 
string later. 

Now, the Dirac strings may end up with monopoles. The action associated with the monopoles is not 
zero at all but rather diverges in ultraviolet, since 

^ J- 

47r e^a 

where a is a (small) spatial cut off. If the length of a closed monopole trajectory is of order L, then the 
suppression of such configuration due to a non- vanishing action is of order 

e"-^ exp {-const L/e^). (71) 

On the other hand, there are different ways to organize a loop of length L. This is the entropy factor. It is 
known to grow exponentially with L as ~ exp ( const' L). 

Thus, one comes to the conclusion (Polyakov, 1975) that at some e = ecrit ~ 1 there is a phase transition 
corresponding to condensation of the monopole loops. As a result, if external electric charges are introduced 
as a probe, their potential energy grows with distance, V(r) ~ r and they are confined. ^ 

To complete the presentation we should explain how one should understand the theory ( |67[ ) that it would 
imply a vanishing action for the Dirac string. 

The crucial point is to define the theory by means of a lattice regularization. Then the action can be 
understood as a sum over plaquette actions: 

S = XI ^ ~ ^^P(* j> A^.dx^') 

X] ^ (1 ~ exp(iF^^,dcr^^,)) = ^-^{^^ cos(F^^dcr^^)) , (72) 

where the sum is taken over all plaquettes. In the continuum limit one reproduces of course the action (|6^). 
However, from the intermediate steps it is clear that the action admits for a large jump in F^^^,: 

F^. ^ F^, + 2TT6i<j^,) (73) 

where the (5- function on the surface is defined as 6{(7fi^)da^'^ — 1 (no summation over /i, i^). The second term 
in Eq. ( |7^ ) exactly corresponds to the Dirac string. Thus, the Dirac strings have no action in the compact 
version of the U{1) gauge theory. 

Note that the UV scale 1/a is the only scale at the model. Thus, it rather exists only as a lattice theory. 
The confining potential sets in for e > Ccrit at all the distances and in this sense the situation does not 
imitate QCD. There is a possibility that a truly continuum theory can be defined near e — Ccrit but we 
would not discuss this topic here (for a recent discussion and further references see, e.g., Jersak et.al (1999)) 



Topology of gauge fixing 

The phenomenon of the monopole condensation is well established in QCD through numerical simulations 
on the lattice, for review see (Chernodub et. al., 1998), (Di Giacomo, 1998).. 

The path to understanding magnetic monopoles in gluodynamics necessarily goes through the paper by 
't Hooft (1980) where the monopoles were introduces as local objects and on general topological grounds. It 
is convenient to start, however, with the Abelian Higgs model, see Eq. (|57|), which contains rather strings 
than monopoles, to learn a new technique. The gauge transformation rotates the phase of the charged field 

$ e"^$. (74) 

3 This confinement mechanism and confinement mechanism of magnetic monopoles in AHM discussed in the previous section 
are interrelated. It can be shown in the lattice regularization that compact U{1) theory is dual to AHM in which the bare Higgs 
and gauge boson mass are infinite. 
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In particular one can rotate in such a way that 



/m$ = (75) 

which defines what is caUed the unitary gauge. The reason for the name is that if < $ >^ then the 
excitation spectrum in this gauge consists of physical massive vector and scalar particles, with no ghosts. 

The crucial point is that, as remarked by 't Hooft, the condition ( [ts]) fixes the gauge uniquely unless 
the real part of $ is also vanishing. Thus, there are exceptional points characterized by the equations 

Ira $ = 

i?e $ = 0. (76) 

Viewed as two equations in four dimensional world, these conditions define a world sheet, or trajectory of a 
string. Note that the definition ( |7^ ) is local in the sense that it defines an infinitely thin string. 

The definition does not give any immediate clue as to whether the strings (|7^) arc important dynamically. 
However, one may notice that <& does vanish on the central axis of the ANO string, see (|6l|). If one makes 
such an identification then it is easy to conclude that the world sheets ( |7^ ) are either closed, corresponding 
to closed ANO strings, or end up with monopoles. 

In case of a nonabelian theory, one introduces first the so called U{1) projection which is nothing else 
but (partial) gauge fixing. Namely, one chooses a field in an adjoint representation. If we concentrate of 
SU{2), then this field is a triplet <I>°. The field $° could be a component of the nonabelian field strength 
tensor F^^ or a composite operator. Then one chooses the gauge to rotate <!> = to, say, the third 

direction in the color space: 

0-3 

$ -> nm-^ " 

where 51 is the matrix of the gauge transformation. The condition ( [77[ ) fixes the gauge only up to rotations 
around the third axis in the color space, that is up to a remaining U{1) symmetry. 

What is crucial for introduction of the monopoles is that the condition ( [tT] ) can be implemented every- 
where except for the points 

= 0. (78) 

Thus, the exceptional points are given by solutions to three equations in the four-dimensional space and 
specify therefore a trajectory. Actually, this is a monopole trajectory. To prove this, one follows exactly 
the same sequence of steps as in establishing the relation of the 't Hooft-Polyakov monopole to the Dirac 
magnetic monopole, see ('t Hooft, 1974). Indeed, in Georgi-Glashow model we diagonalize the higgs triplet 
<I>° as in ( [ttI) and j'l'p vanishes at the center of a 't Hooft-Polyakov monopole. And the basic observation is 
that the algebra associated with the rotation of the 't Hooft-Polyakov monopole to the Dirac monopole is in 
fact uniquely determined by this circumstance. A detailed derivation can be found, e.g., in (Simonov, 1996) 
or in (Chernodub et.al., 1998). 

Note that the definition of the monopole (^8|) is pure topological, the same as the definition of the string 
above, see (|76|). It does not at all imply existence of a classical monopole solution. However, it provides with 
a local definition of a magnetically charged field ^m{x) and allows to study the monopole condensation. 
Note also that the definition of the monopoles depend on the choice of the auxiliary field (operator) <I>° . The 
dual superconductor model of confinement is confirmed numerically most impressively within the so called 
maximal Abelian gauge. In this case, one fixes gauge, up to a U{1) subgroup by minimizing the value of 
(A^)^ -I- (A^)^, where Aj^'^ are components of the vector potential in the color space. The remaining U{1) 
symmetry is the phase rotation of the complex field Aj^ + lA^^. We note all this only in passing, and the 
reader is advised to consult (Chernodub et. al., 1998) for any detail. 

To summarize: both strings in the Abelian Higgs model and monopoles in gluodynamics can be defined 
locally in topological terms. These definitions do not provide however, with any straightforward way to 
evaluate the dynamical significance of these objects. 
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Ultraviolet regularization and non-perturbative effects. 

Now we would like to introduce a new idea on the connection of the topological strings and monopoles 
discussed in the preceding section and nonperturbative ultraviolet divergences. Indeed, intuitively it is 
appealing to guess that, if the string- or monopole-like solutions with finite energy exist only for theories 
with scalar fields, similar point-like, or topological excitations build up on the gauge fields alone have infinite 
action, generally speaking. Only if we accept a prescription of an ultraviolet regularization which eliminates 
these divergences, such excitations may play a dynamical role. Although intuitively the idea looks appealing, 
we shall not be able to prove it in its generality but rather illustrate it on examples. 
Let us begin with the Dirac string. Naively, the energy eo is equal to: 



U'd'r ^ ^-^^^^A I dl, (79) 



where A is the cross section of the Dirac string. Thus, eo diverges as 1/A (quadratically) in ultraviolet once 
A^O. 

Instead, we assumed ££> = 0. There are two ways to justify this assumption: 

(i) Lattice regularization, as discussed in the preceding subsection. 

(ii) Duality between electric and magnetic charges. Since the electric charges do not have strings at- 
tached, the action associated with the Dirac string is postulated to vanish. 

Note that while in case of the compact U{1) we relied on point (i), while in case of the Abelian Higgs 
model we had to rely on (ii). 

Turning to SU{2), there is of course a possibility to introduce a Dirac string associated with any U{1) 
subgroup. In this connection, let us mention another language for the strings, that is singular gauge trans- 
formations, for details see, e.g., Chernodub (1998). Gauge transformations can be specified in terms of a 
matrix fi: 

n{x) = exp ^ia°(2:)y^ (80) 
where a°'(a — 1, 2, 3) are parameters of the transformation and cr" are the Pauli matrices. Then 

n^{x) (^i^(x) - ^d^^ n{x) (81) 

where = yl°(T''/2. Moreover for the field strength one has 

F,,, ^ F^M - ^n\d^d, ~ d,d^)n (82) 

where we reserved for the possibility that the derivatives 9^,,^ do not commute when applied to fl{x). The 
noncommutativity is in fact the definition of a singular gauge transformation. 
Consider furthermore 

- ( ''i.. ""^O (83) 
— sm^e cos^ J 

where a and 7 are azimuthal and polar angles, respectively. Then it is straightforward to check that 

27r 

F^,{A'') = n^F^,{A)n - — ('5^,1^-2 - 5^^2Su,i) Sixi)Six2) e{-xs) (84) 

In other words, we generated a Dirac string directed along the xs-axis ending at 0:3 = and carrying the 
color index a = 3. 

It is quite obvious that such Abelian-like strings are allowed by the lattice regularization of the theory. 
However, we cannot get too far using only Abelian-like field configurations. Indeed, imagine that the string 
(Q) would go over into the same Abelian monopoles at its end points. Then we would have the same estimates 
as in the case of the compact U{1) but our coupling now tends to zero in the ultraviolet, \ima^ag^{a) = 
and these field configurations are strongly suppressed. 
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The best qualitative picture for QCD based on the analogy with the compact U{1) seems to be as follows. 
Begin with a small lattice size a but then go to a coarser lattice, with corresponding renormalization of the 
coupling a la Wilson. Then once the effective coupling governing a U{1) subgroup reaches the value of e^^jj 
then one may think that the monopoles are condensed, in the same way as in the U{1) case. There is a 
semi-quantitative prediction based on this picture. Namely, it is natural to assume that the running of the 
coupling is changed drastically at this point as far as further advance to the infrared is concerned. Thus, we 
can speculate that the coupling is frozen at the value 

5/rozen ~ '^^crit ~ ^ (85) 

which is not unreasonable phenomenologically^. 

To summarize, we touched upon a few issues in this subsection. First of all, we emphasized that 
introduction of the Dirac string implies, at least naively, a new quadratic ultraviolet divergence. We also 
argued that the Abelian-like monopoles cannot be important in the non-Abelian case since the bare coupling 
goes to zero. In a way, the non-Abelian monopoles are to be so to say empty at short distances in the sense 
that the action is determined by contribution of distances of order Aqcd^^ where the coupling is not small. 
How this might happen, we shall explain in the next section. 



Magnetic monopoles in the limit mn oo. 

Magnetic monopoles of finite energy are found ('t Hooft (1974), Polyakov (1974)) as classical solutions to 
nonabelian field theories which include scalar fields. Thus, it is a common question, how one can understand 
monopoles in QCD when there are no scalar fields. We will provide a partial answer to this question, arguing 
that the scalar fields are replaced in some sense by singular gauge transformations as far as short distances 
are concerned. 

Begin with the Georgi-Glashow model which includes a tripled of scalar fields: 

L = liPp' + liD,<i>'^r + Inm' - (86) 



Here, 



= d,Al-d,Al + ge'^>'^AlA: (87) 
Df,^" = a^*'' ge^^'^Aj^^^ (88) 

and A, v are constants. 

Furthermore, one introduces functions A{r), <i>(r): 

j^a ^ £fab/^(r), = r''$(r). (89) 

The monopole solution is characterized by its asymptotic at r — > oo: 

\im A{r)^^^ = -r-^g-\ lim$(r)^^^ = v r'K (90) 

The mass of the monopole is: 

M^onopoie = C{X/g^), (91) 

9 

where C{X/g'^) is in fact a very slow function of its argument, as emphasized already by 't Hooft (1974). 
Eq. (|9l] ) implies then that the mass of the monopole corresponds to the mass of a Dirac monopole with an 
ultraviolet cut off at distances of order r ~ rriy^ . 

What is a crucial point for us now, is that the mass of the monopole remains finite in the limit A — > oo 
which is the same as the limit of infinite Higgs mass, rnj^ = \v. Thus, we may take the scalar mass as an 

''Note that the lattice calculations show (Kato 1998, Chernodub 1999) that the abelian monopoles of finite physical size 
("extended" monopoles (Ivanenko 1990)) are important for confinement. 
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ultraviolet cut off and hope to remove the scalar from the spectrum altogether. Now we are dealing with 
QCD still having a finite mass of monopole. Let us look closer (in the limit mn 3> my) why we do need the 
scalar field. 

Consider distances niy^ 3> r 3> The point is that at such distances the field of the monopole can 

be nothing else but nearly a pure gauge, 

« 0. (92) 

Indeed, only in this case the mass of the monopole would be insensitive to distances much smaller than rriy^ , 
as testified by (|9l|). To find the solution more explicitly, let us look closer at the equations of motion. Since 
we consider r ^ "i^^^ we choose $(r) in its asymptotical form, r — > oo, see (|90|). Moreover, we make an 
ansatz: 



c 



A{r) « — (93) 
where c is a constant. 

Then the equations of motion reduce to a simple algebraic equation (see, e.g., 't Hooft (1974)): 

= Acg-^ + 6g-^c^ + 2g-^c^ + 2v^r^gil + c). (94) 

At the distances we are considering the product gvr <^ 1 and let us neglect for the moment the last term in 
(H). Then we have three solutions for the constant c: 

ci = -1, C2 = -2, C3 = 0. (95) 

The first solution is singled out by the fact that it nullifies also the terms in (|9^) . This is crucial at large 
r and that is why ci is valid asymptotically, see (|90|). However, this solution implies (Fl^^)'^ ^ r^^ and a 
large contribution to the mass from distances r <^ rriy^ which is not allowed, see (^l]). On the other hand, 
if we choose C2 = — 2 then Eq. (^2|) is fulfilled and the contribution to the action from the non-Abelian field 
vanishes. Thus, at small distances we should turn rather to the solutions 02,3. To remove the contradiction 



with the Eq. ( |94| ) we include the next term in the expansion (93) 



A{r) w JL+dv^lnr, (96) 

gj,Z 

where d is a constant and the second term is a small corrections at the distances we are considering. The 
constant d can be readily found: d = ± g/3. 
Since ~ 0, the potential 

is obtainable from A° = by a gauge transformation: 

- 4^''°'^'? = i{^ydf,il°. (98) 

It is not difficult to find out the explicit form of the matrix 

Vl^ = ia-n (99) 

where n is the unit vector from the center of the monopole to the observation point. 

Now, it is clear that if we would introduce gauge rotations (p9| ) ad hoc, without scalar fields, then we 
would conclude that the corresponding A^ has a (5-function singularity at the origin so that not only 
but the total mass ~ / (F^y)^d"^r is UV divergent. In this sense the only function of the scalar field is to 
smoothen the fields at the origin so that the contribution of the singularity is in fact normalized to zero. 
Now, the observation is that if we use the lattice regularization and ask the question whether the singular 
(at r = 0) potentials ( ^ ) are allowed or not, the answer would be in positive. Namely, the center of the 
singular monopole would fall into a center of a lattice cube and this would imply that the singularity does 
not contribute to the action in fact. 
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Short Strings. 

We proceed now to derivation of the stringy correction (|56| ) to the potential V{r) within the Abehan Higgs 
model (Gubarev et. al., 1998). Thus, we consider the following problem. Two external magnetic charges are 
put at a distance r in the vacuum of the Abelian Higgs model. In particular, the vacuum expectation value 
of the scalar field <f> is non- vanishing. The problem is to find the potential energy V{r) at distances much 
smaller than mjj^ , nriy^ . 

The energy is determined in the classical approximation and, at first sight, the problem is not much of a 
challenge. The crucial point is that one has to impose an extra boundary condition, namely vanishing of the 
scalar field along a mathematically thin line connecting the magnetic charges. We have already mentioned 
this condition while discussing the topological strings (see Eq. and the related discussion). 

The language of the Dirac string can be useful to substantiate the point. Indeed, the infinitely thin 
line discussed above is nothing else but the Dirac string connecting the monopoles. The possibility of its 
dynamical manifestations arises from the fact that the Dirac string cannot coexist with <i> 7^ and <i> vanishes 
along the string. Indeed, self-energy of the Dirac string is normalized to be zero in the perturbative vacuum. 
To justify this one can invoke duality and ask for equality of self-energies of electric and magnetic charges. 
Since the electric charge has no string attached, the requirement would imply vanishing energy for the Dirac 
string. However, if the Dirac string would be embedded into a vacuum with ($) ^ then its energy would 
jump to infinity since there is the term l/2|$pA^ in the action and 00 for a Dirac string. Hence, 

$ = along the string and it is just the condition mentioned above. In other words, Dirac strings always rest 
on the perturbative vacuum which is defined as the vacuum state obeying the duality principle. Therefore, 
even in the limit r — *■ there is a deep well in the profile of the Higgs field <&. This might cost energy which 
is linear with r even at small r. 

The next question is whether this mathematically thin string realizes as a short physical string. Where 
by the physical string we understand a stringy piece in the potential, (Tor at small r. In other words, we are 
going to see whether the stringy boundary condition implies a stringy potential. To get the answer we solve 
classical equations of motion. Let us note that the scheme of the calculation and some numerical results 
for the potential V{r) can be found in a number of papers, see, e.g., (Ball and Caticha, 1988). However, 
prior to (Gubarev et. al., 1998) there were no measurements dedicated specifically to small corrections to 
the Coulombic potential at r ^ 0. 

We will consider the unitary gauge, /m$ ~ 0. Then the most general ansatz for the fields consistent 
with the symmetries of the problem is: 

$ = rif{p, z) Imf = 



(100) 



1/2 



Z — X3 Xa — -^aj P 



a 1,2 



In the limit r — > the Coulombic contribution becomes singular. The easiest way to separate the singular 
piece is to change the variables A = Ad + a, where Ad is the solution in the absence of the Higgs field: 



Ad = 



1 



£+ 



z± = z ± r/2 r± = [p^ 



1 1/2 



(101) 



Let us introduce also a new variable k = \/2X/e = mu /my and measure all dimensional quantities in terms 
of my — erj. Upon the change of variables ( |lOl| ) the energy functional and the classical equation of motion 
take the form: 



E{r) = ^ J dz J pdp\[^dp{pa)]^ + [dM'+ 
-00 ^ 

+ [dpf? + [dj? + P{a + Ad? + \n\p - if] 



(102) 
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'-dp{pa) 



+ dla = f{a + Ad) 



-dp [pdpf] + dlf = f{a + + l-K^fif - 1) 



(103) 



(104) 



The energy functional has been minimized numerically. The numerical results (Gubarev ct. al., 1998) 
for various k values clearly demonstrate that there is a linear piece in the potential even in the limit r <C 1. 
The slope ctq at r ^ was defined by the fitting the numerical data to: 



Efitir) = Co 



1 



-1 ) + K + ^Co)r = <7or + 0(r2), 



(105) 



The resulting slope ctq depends smoothly on the value of k, For the purpose of orientation let us note that 
for K = 1 the slope of the potential at r — > is the same as at r — > oo That is, within error bars: 



(106) 



where a^o determines the value of the potential at large r. 

To summarize, existence of short strings has been proven in the classical approximation to the Abelian 
Higgs model. The linear piece in the potential at small distances reflects the boundary condition that $ = 
along the straight line connecting the monopoles. 



Strings vs OPE. 

Let us emphasize that both in case of the compact U{1) and in case of the Abelian Higgs model the OPE 
does not apply to study the power corrections. 

In case of the compact U{1) it is especially evident that the OPE does not work. Indeed, as we emphasized 
many times, the use of the OPE allows to parameterize infrared sensitive parts of the Feynman graphs. But 
there are no Feynman graphs at all now because we are considering the Lagrangian of a free gauge field. 
Nevertheless, the answer depends on the value of the coupling e^. If is higher than the critical value 
e^^jj ~ 1 then the system is in confinement phase (see above). As a result the physics is changed at all the 
distances. Indeed, the ultraviolet scale a (which may be thought of as the lattice size) is the only scale in the 
problem since the coupling does not run. As a result, the potential between external electric charges becomes 
linear at all the distances. One photon exchange is not seen at all. Clearly, this observation is equivalent to 
saying that the OPE does not work at any distance because the OPE is just the assumption that at short 
distances one sees particle exchange. Thus, we may say that the example of the compact U{1) is too strong. 
If > ej?,.jj then there is no particle exchange at all. While in QCD we believe that perturbation theory 
applies to describe the leading effects at large Q^. 

Turn now to the AHM where we can check OPE against the direct calculation of the potential in the 
preceding subsection. At first sight, existence of the l/Q^ corrections in case of the AHM is not surprising 
since now there is an operator of dimension d = 2, that is It is worth emphasizing therefore that a bit 
more careful analysis demonstrates that the power correction to the potential which we found in the AHM 
also contradic;ts the OPE. 

Moreover, it is the existence of short strings that is manifested also through breaking of the standard op- 
erator expansion. Indeed, above we found the potential in the classical approximation. In this approximation 
the potential is usually directly related to the propagator Di^v{(f) in the momentum space, 

V{r) = / e^'*" ^oo(q') (107) 

Moreover, as far as is in the Euclidean region and much larger than the mass parameters, the propagator 
D^u{q^) can be evaluated by using the OPE. Restriction to the classical approximation implies that loop 
contributions are not included. However, vacuum fields which are soft on the scale of can be consistently 
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accounted for in this way (for a review see, e.g., Novikov et.al. (1985)). This standard logic can be iUustrated 
by an example of the photon propagator connecting two electric currents. Modulus longitudinal terms, we 
have: 

D,M') = <5,.fl + + + - (108) 
\q q q q q q J q — my 

Thus, one uses first the general OPE assuming \q^ \ ^ e^^^ then substitutes the vacuum expectation of the 
Higgs field $ and upon summation of the whole series of the power corrections reproduces the propagator 
of a massive particle. The latter can also be obtained by solving directly the classical equations of motion. 

This approach fails, however, if there are both magnetic and electric charges present. In this case, one 
can choose the Zwanziger formalism (Zwanziger, (1971), Brandt (1978)) and work out an expression for 



propagation of a photon coupled to magnetic currents following literally the same steps as in (108). In the 
gauge n^D^u ~ the result is well known: 

q'^ — niy \ (qn) ^"^^ ■ -if /j/ ' (^^^^2 ' (g/j)^ 



Dp.i,iq,n) = — — — ( ^^.f - TTTY^^A'"-'^ + gi^"-^) + /!!„\2 + f„J\2 i^t^'^'^^ ~ "p"'^) ) ■ (109) 



Here the vector is directed along the Dirac strings attached to the magnetic charges and there are general 
arguments that there should be no dependence of physical effects on n^. On the other hand, if the potential 



energy is given by the Fourier transform of (109) then its dependence on is explicit. 

Note that Eq. (109) immediately implies that the standard OPE does not work any longer on the level 
of corrections. Indeed, choosing q^ large and negative does not guarantee now that the niy correction 
is small since the factor (qn)^ in the denominator may become zero. Of course, appearance of the poles in 
{qn) in the longitudinal terms is not dangerous since they drop due to the current conservation. However, 



the term proportional to my in (109) cannot be disregarded and contribute, in particular, to (107) 



The reason for the breaking of the standard OPE is that even at short distances the dynamics of the 
short strings should be accounted for explicitly. In particular, in the classical approximation the string lies 
along the straight line connecting the magnetic charges and affects the solution through the corresponding 
boundary condition, see above. More generally, the OPE allows to account for effects of vacuum fields, in 
our case for ($) ^ 0. The OPE is valid therefore as far as the probe particles do not change the vacuum 
fields drastically and the unperturbed vacuum fields are a reasonable zero-order approximation. In our case, 
however, the Higgs field is brought down to zero along the string and this is a nonperturbative effect. Thus, 
the stringy piece in the potential V{r) at r — > is a nonperturbative correction which is associated with 
short distances and emerges already on the classical level. 



4 REVISITING THE PHENOMENOLOGY. 
Revisiting the quark potential at short distances. 

We have shown above that in the AHM the quark potential contains a linear correction at short distances. 
Moreover, the reason for the breaking of the OPE seems to be a unified description of magnetic and electric 
charges. Since the monopole condensation is established within the abelian projection of QCD, we would 
speculate that a similar picture for the potential is true in QCD as well. Then we expect a linear piece 
dominate over other non-perturbative corrections at short distances and this would make it unnecessary to 
account for the retardation effects discussed in part I. 

Because the theoretical foundation for the new picture is not absolutely solid, phenomenological evidence 
would be especially important at this point. 

Without going into detail, because of the space consideration, let us mention that there are a few areas 
where measurements, especially on the lattice, could clarify the situation: 

(«) Direct measurements of the potential V{r) bring so far nonvanishing results for the effective string 
tension at short distances (Bali et al, 1995, Bali, 1999): 

(To = (1^6)aoo (110) 
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{ii) There exist lattice measurements of fine splitting of the QQ levels as function of the heavy quark 
mass |. The Voloshin-Leutwyler picture results in a particular pattern of the mass dependence of this 
splitting. Moreover, these predictions are very different from the predictions based on adding a linear 
potential to the Coulomb potential at short distances. Numerically, predictions based on short strings 
with ctq w CToo (see above) are close to the predictions obtained by BuchmuUer and Tye (1980) within a 
phenomenological model for the potential. 

The results of most advances measurements of this type are published in (Fingberg, 1998) and favor 
strongly the linear correction to the potential at short distances. 

(iii) There is an interesting evidence that the nonperturbative fluctuations on the lattice responsible for 
the confinement can be identified as the so called P vortices (Faber et. al., 1999). If one measures the quark 
potential due only to these vortices, the numerical results seem to indicate that the slope of the potential is 
the same at large and small distances, see (L. Del Debbio et. al., 1997). 

(iv) Analytical studies of the Bethe-Salpeter equation and comparison of the results with the data about 
the charmonium spectrum favor a nonvanishing linear correction to the potential at short distances (Badalian 
and Morgunov, 1999). 

Thus, it appears that existing experimental data favor a linear correction to the heavy quark potential 
at short distances but much more work is needed to finalize the results. 

To summarize: the stringy correction to the Coulombic potential at short distances found in the Abelian 
Higgs model violates the operator product expansion. The search for a linear piece in the potential at short 
distances in QCD is of great importance. We listed above some preliminary indications that such a correction 
is indeed present. 



Revisiting correlation functions: tachyonic gluon mass at short distances. 

The l/Q^ corrections discussed in the previous section go beyond the standard OPE. Detection of the 
new type corrections through phenomenological studies would be of great interest. In this section we will 
discuss phenomenology in terms of a tachyonic gluon mass which is assumed to mimic the short-distance 
nonperturbative effects (Chetyrkin et al., 1998). 

First, let us note that not all the l/Q^ corrections in QCD are associated with short distances. For 
example, in case of DIS the corrections are coming from the IR region and perfectly consistent with 

the OPE. Thus, the class of theoretical objects for which an observation of the corrections would signify 
going beyond the OPE is limited. One example was the potential V{r) discussed above. Other examples are 
the correlator functions (|^) where the IR-sensitive power corrections start with terms, as is explained 
in detail in the Part I of the lectures. 

Thus, we concentrate on this set of variables and an interesting question is whether there is room for 
introduction of sizable non-standard corrections to the correlator functions. The answer seems to be 

in positive. 

Even if one accepts that the non-standard correction to the potential has already been observed 

(see the preceding section), it is far from being trivial to evaluate the corrections to other quantities. 

Qualitatively, however, one may hope that introduction of a tachyonic gluon mass at short distances would 
imitate the effect of the ^qcd/Q^ corrections. Indeed, the linear term aQ-r in the potential at short distances 
can be imitated by the Yukawa potential with a gluon mass A: 

- - ao. (Ill) 

6 



On the Born level, Eq. (Ill) is an identity. However, a tachyonic gluon mass can be consistently used at 



one-loop level as well. Of course Eq. (Ill) may serve only for a rough estimate 



Looking for other applications, let us remind the reader, that one of the basic quantities to be determined 
from the theory is the scale M^^^j at which the parton model for the correlators gets violated considerably 
via the power corrections. Now, a new term proportional to A^ is added to the theoretical side of Ilj{M^) 



^ We are thankful to A. Leonidov for bringing this issue to our attention and for introducing to the literature. 
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which becomes: 



(112) 



where Cj is calculable in terms of (Chetyrkin et.al., 1998)^: 

4q!, 4a 



4Cp — Cqluonium A . (113) 

OTT TT 



Phenomenologically, in the p-channel there are severe restrictions (Narison, 1992) on the new term Cp/M^: 

Cp sa - (0.03 ~ 0.07) GeV^ (114) 



Remarkably enough, the sign of Cp does correspond to a tachyonic gluon mass (if we interpret Cp this way) 



p 

Moreover, when interpreted in terms of the constraint (114) does aUow for a large A^, say, A^ = —0.5 GeV 



As for the 7r-channel one finds now a new value of M^^^j associated with A^ 7^ 0: 

^criti'^ - channel) « 4 • M^^^^{p - channel). (115) 

It is amusing that just this value of M^^^^ was found in (Novikov et.al., 1981) from an analysis of the pion 
contribution and it cannot be explained within the standard QCD sum rules. Moreover, the sign of the 
correction in the 7r-channel is what is needed for phenomenology (Novikov et.al., 1981). Fixing the value of 
Ctt to bring the theoretical n7r(M^) into agreement with the phenomenological input one gets 

A^ « -0.5 GeV^ (116) 

Finally, with this value of A^ in hand, we can determine the new value of A/^^j^ in the scalar-gluonium 
channel and, again, it turns to be what is needed for the phenomenology, see Eq. (^). 

Further checks could be provided by measurements of various correlation functions on the lattice, similar 
to the measurements reported in (Chu et. al., 1993). It would be most interesting to try to disentangle the 
corrections discussed here from the effects of the direct instantons mentioned in Part 1. Let us give 
two particular examples: 

(z) In case of instantons, the correlation function of the tensor currents, which essentially coincide with 
the energy-momentum tensor are protected against large corrections. Indeed, it is well known that the energy- 
momentum tensor vanishes on the instantons. If the tachyonic mass is responsible for large corrections, the 
tensor gluonic current is not singled out in any way. 

(ii) It would be very important to measure the sum 

Il„{x)+n^{x) 

at critical values of x where the violation of the asymptotical freedom becomes noticeable. The point is that 
the direct instantons drop off from this sum to the first order while the effect of the tachyonic mass would 
add up. The existing data (Chu et. al., 1993), taken literally leave ample space for the tachyonic mass. 
Measurements with improved accuracy are highly desirable. 

It is worth emphasizing also that the A^ terms represent nonperturbative physics and limit in this sense 
the range of applicability of pure perturbative calculations. This nonperturbative piece may well be much 
larger than some of the perturbative corrections which are calculable and calculated nowadays. 

To summarize: at least qualitatively, the phenomenology with a tachyonic gluon mass which is quite 
large numerically stands well to a few highly nontrivial tests. Further crucial tests of the model with the 
tachyonic gluon mass could be furnished with measurements of various correlation functions Ilj{]VP) on the 
lattice. 



Further discussion of the short-distance tachyonic gluon mass can be found in (Zakharov, 1998), (Simonov, 1999), (Huber 
et.ah, 1999). 
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Revisiting vacuum condensates: condensates of dimension = 2 in gauge theories. 



As is explained in the first part of the lectures, introduction of various vacuum condensates turned a useful 
way to understand and characterize dynamics of QCD. The most famous example is the quark condensate: 

{0\qq\0) ^ (117) 

where q stands for light u- or d-quarks. In the realistic case of negligible small quark mass, a nonvanishing 
value of the quark condensate manifests spontaneous breaking of the chiral symmetry. 
Within the framework of QCD sum rules there emerged also the gluon condensate 

(0|a,(G«J»^0 

A non- vanishing value of this condensate docs not signify breaking of any symmetry but rather provides with 
a measure of nonpcrturbative fields in the vacuum. The gluon condensate appears to be the simplest gauge 
invariant quantity characterizing nonpcrturbative vacuum fields. It has dimension d = A and this dimension 
is echoed in the observation that the leading non-perturbative corrections in the QCD sum rules at large Q 
are of order (0|q;s(G°^)2|0)/Q*. 

Now, wc are driven to a rather astonishing conclusion that condensates of dimension d = 2 may also 
have dynamical significance. Indeed, the condensation of monopoles in QCD is described by < <1?^ >^ 
which is a dimension d = 2 vacuum expectation value. Also, we saw that in the compact U (1) there arises a 
nonvanishing string tension. The condensates are of nonpcrturbative nature and their appearance is related 
to the compactness of the gauge group. Clearly, such condensates cannot be detected by means of the 
standard OPE since they are not related to softness of propagators in the Feynman graphs. 

There is no direct clash with the gauge invariancc since <1>m does not enter the Lagrangian but is 
constructed in a rather indirect way, see the preceding section. The situation when there emerge condensates 
not related at all to the original fields may appear, however, confusing. In terms of the original fields 
condensates of dimension d = 2 can be related only to the vacuum expectation value of < Aj^ >. And 
in this section we will argue, following (Gubarev et. al., 1999) that the condensate < > constructed 
directly on the vector-potential of the gauge field may also have dynamical significance, upon some 
clarifications Literally such a condensate is gauge non-invariant and devoid therefore of physical meaning. 
We will consider, however, the minimal value of it, < {A^)min >• Looking for the minimal value of A"^ 
implies, of course, a particular choice of the gauge. Also, we have to consider the Euclidean space-time so 
that the minimal value of could be meaningful. 

The relevance of < {A^)min > can be understood on a simple example of very thin cosmic strings 
borrowed from (Alford and Wilczek, 1989). Then the field strength tensor cannot be detected directly. 
However, because of of the Aharonov-Bohm effect particles can be scattered off the strings. From our 
point of view, it is crucial that the system is characterized by a non- vanishing A^. Indeed, because of the 
topological condition, 

^A|,dx^' = j H-rfs = (118) 

where is the magnetic flux, the value of A^^ cannot be brought down to zero everywhere. Thus, the value 
of < {A^)min distinguishes the system from the trivial vacuum for which, on the classical level, one 
can choose = 0. 

In case of the compact U{1) gauge group, one can measure the < {A'^)min > numerically. Moreover, it 
is possible to substruct explicitly the perturbative contribution to this quantity leaving only the dynamically 
relevant part < {Afjrnin >non-pert.- The result for < {Af^)min >non-pert. is presented in Fig. 6. We see 
that the behavior of this condensate clearly signals the phase transition to the confining phase near = 1. 
Although, it cannot be considered, strictly speaking, as an order parameter since even the nonpcrturbative 
part of < (A^)toto > does not vanish in any phase. 

To summarize: dimension two vacuum expectation values, like < {A^)min > may have a hidden topo- 
logical meaning and be dynamically significant. 
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Figure 6: Non-porturbative contribution to < > as a function of /3 = 1/e^. The plot is borrowed from (Gubarev et. 

al., 1999) 



5 CONCLUSIONS. 

In Part I we considered power corrections within the OPE approach, which culminates in the QCD sum 
rules (see, e.g., (^S])). The sum rules allowed to tune the power corrections to resonances in some channels, 
like the p-meson channel. There are many further successful applications, see, e.g., Reinders et.al. (1985), 
Narison (1989), which we just did not have time at all to cover. More recently the technique was generalized 
to observables directly in the Minkowski space by using infrared renormalons or infinitesimal gluon mass. 

On the other hand, in some so to say exceptional channels the sum rules certainly fail (Novikov et.al., 
1981). Most conspicuous cases are the 7r-meson, 0"'"-gluonium channels. There are some explanations to 
the failures in terms of direct transitions of the currents to instantons (Novikov et.al. (1981), Shuryak and 
Schafer (1998)). However, the instanton corrections are lacking simple analytical structure. Also, neither 
the soft vacuum fields accounted for by the OPE nor the direct instantons seem to encode the effects of 
confinement. 

Moreover, very recent studies of the power corrections in the Abelian Higgs model reveal novel 1/Q^ 
corrections which are associated with small distances and are not accounted for in the operator product 
expansion. If such corrections exist in the QCD case as well, then there are quite a few places where the 
new corrections could be detected through measurements, mostly through measurements on the lattice. The 
existing data are not discouraging for the new corrections, not at all. Further checks are desirable. 

To summarize, the nature may turn to be generous as far as power corrections are concerned. I am 
borrowing this term from a talk on dark matter. Indeed, first people assumed that there should be a single 
dominant source of the dark matter, and now it appears distributed among various equally important com- 
ponents. Similar picture may be true for the sum rules. Indeed, very first idea would be that theoretically 
the correlation functions IIj{Q'^) could be found perturbatively at large and the growth of the effective 
coupling at smaller would signal the breaking of the asymptotic freedom. Then the picture got more 
involved and the effect of soft non-perturbative fields was included in terms of the quark and gluon conden- 
sates. It appears not suffice to explain the peculiarities of all the channels and the effect of direct instantons 
was invoked. As the latest development, hypothetical l/Q^ corrections associated with short strings are 
established within the Abelian Higgs model which is thought to mimic the QCD confinement. 

Thus it appears now that practically all known "ingredients" of QCD have already found their way into 
the physics of the power corrections. 
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